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Abstract 

Through the scattering of light-pseudoscalar mesons (n, K, r), r}') off charmed mesons 
(D,D S ), we study the D* 0 (2317) state and other relevant charm scalar mesons in a unita- 
rized chiral effective held theory approach. We investigate the charm scalar meson poles 
with different strangeness (S) and isospin (J) quantum numbers as well as their correspond¬ 
ing residues, which provide the coupling strengths of the charm scalar mesons. Both the 
light-quark mass and Nc dependences of the pole positions of the D* 0 (2317) and the poles 
with ( S, I) = (0,1/2) are analyzed in detail in this work. Interestingly we observe quite sim¬ 
ilar pion mass trajectories for the resonance pole at around 2.1 GeV with (S , I) = (0,1/2) to 
those of the /o(500) given in the literature. When increasing the values of Nc we find that 
a bound state and a virtual state in the ( S,I ) = (1,0) channel asymmetrically approach 
the DK threshold for Nc < 6, and they meet at this threshold at Nc = 6. When Nc > 6, 
the bound and virtual states move into the complex plane on the second Riemann sheet 
and become a symmetric pair of resonance poles. For large enough values of Nc, neither 
the D* 0 (2317) pole nor the poles with (S, I) = (0,1/2) tend to fall down to the real axis, 
indicating that they do not behave like a standard quark-antiquark meson at large Nc- 
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1 Introduction 


The discovery of the scalar charmed-strange meson D* 0 (2317) [TH3] has triggered intensive 
studies in hadron physics. With its mass just about 40 MeV below the DK threshold, it is 
one of the few precious examples of a loosely bound state observed in the meson sector. The 
scattering of the charmed and light pseudoscalar mesons provides a powerful tool to probe the 
inner structure of D* 0 (2317). Recently much progress from lattice simulations focusing on the 
D* 0 (2317) state has been made investigating the previously mentioned scattering processes [4- 
j7j. The pion masses used in these lattice simulations are still large compared to the physical 
value, so the chiral extrapolations are necessary to obtain physical predictions from the lattice 
data. Chiral perturbation theory (yPT), as an effective field theory of low energy QCD, has 
been demonstrated to be a successful and reliable tool to perform the chiral extrapolations in 
many lattice simulations [8j. 

In traditional heavy meson xPT (HMyPT), proposed by Refs [51411]. not only is the chi¬ 
ral symmetry encoded but also the heavy-quark symmetry is complemented so as to study 
heavy-light systems. Due to the incorporation of heavy-quark symmetry, the interactions are 
independent of heavy-quark spin and flavor, which imply that the scalar D and vector D* 
should be treated as counterparts to form, for instance, a spin doublet. However, following 
Ref. [12] , a covariant HMyPT, in combination only with the chiral symmetry, will be imposed 
to investigate the interactions between the charmed D and pseudo Nambu-Goldstone bosons 
(pNGBs) in this work. In covariant yPT, the necessity of the inclusion of the D* relies on 
the size of its contributions. In our case, only the 0 + charmed states involved in the 5-wave 
interactions will be studied and the vector D* only contributes via Ti-channel exchange, which 
provides a negligible contribution to the interacting potential at tree level m- Therefore we 
consider it is a good approximation to exclude the vector D* mesons when focusing only on 
the 5-wave charmed D meson and pNGB scattering. Within the framework of covariant heavy 
meson xPT by taking the and the octet of pNGBs as the dynamical fields, the leading or¬ 
der (LO) [TiUTB] . next-to-leading order (NLO) p~2ll 13 LI 1 7I4T9] and next-to-next-to-leading order 
(NNLO) [201I2T] calculations have been performed to investigate the properties of D* 0 ( 2317). 
The pole trajectories with varying up/down- and strange-quark masses are useful quantities 
to further probe the inner structure of the hadrons, e.g. the M n and Mk dependences of 
D* 0 (2317) have been extensively studied in previous Refs [1311181119] and similar trajectories for 
the light-flavor resonances /o(500)/<t and p(770) are predicted in Refs. [22][23]. Another kind 
of pole trajectories with varying Nc, with Nq the number of colors in QCD, has been demon¬ 
strated to be useful to discriminate different explanations for light-flavor resonances [24H28] , 
such as for /o(500), p(770), etc. However, the study of the Nc trajectories of the _D* 0 (2317) 
state is still lacking. One of the key novelties of the present work is to fill this gap. The 
Nc trajectories for other pertinent charm scalar mesons resulting from the and pNGBs 
scattering shall be investigated as well. We note that finite-volume corrections for this type of 
scattering processes have recently been worked out, see Refs. [3T1I32] . 

We point out that in order to study the Nc trajectories of the poles from the and 
pNGBs scattering, it is inappropriate to only consider the 517(3) octet of pNGBs, i.e. tv, K, rj$. 
The singlet 770 becomes a relevant degree of freedom (d.o.f) when discussing large Nc- The 
reason is that the QCD U(1)a anomaly, which is responsible for the massive % in the physical 
case at Nc = 3, is 1 /Nc suppressed in the large Nc limit. As a result, the QCD spectrum 
is the pNGB nonet at low energy in the chiral and large Nc limits [334436] . rather than the 
517(3) octet resulting in the chiral limit. Therefore it is necessary for us to generalize the 
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previous discussions on the charmed mesons and the pNGB octet scattering, such as 
those in Refs. nsHigimiiaED, to the processes involving the singlet t]q. After the inclusion 
of the singlet '%, not only the QCD spectrum is completed in the large Nq limit, but also we 
can have a more realistic description for the physical meson at Nq = 3, since the physical p 
and rf mesons are mixtures of rjs and r]Q. While in the standard SU(3) yPT, the r)% is typically 
identified with the physical r/ state [37] and rj — rf mixing effects appear in certain LECs, in 
particular in Lj. Within the yPT framework with explicit light pseudoscalar octet plus singlet 
ir,K,r],r]' (recognized as large Nq or {7(3) yPT [38]) and the charmed mesons D^, we give 
an updated discussion by taking into account the lattice simulations and pay special attention 
to the pole spectra in the scattering of light pseudoscalar and charmed mesons. 

The paper is organized as follows. We first calculate all of the amplitudes for the charmed 
meson ( D , D s ) and light pseudoscalars ( 7 r, LT, 77 , 7 /) scattering processes in U(3) yRT, and then 
we use a simple recipe to unitarize the perturbative scattering amplitudes. This will be the 
main focus of Sect. [ 2 ] The fits to the lattice simulation data will be presented in Sect. [2 where 
we determine all of the free parameters. The pole contents, the corresponding residues and the 
pole trajectories with varying light-quark masses and Nq will be analyzed in detail in Sect. [H 
Finally, we give a summary and conclude in Sect. [5] 

2 Chiral amplitudes and their unitarization 

Due to inclusion of the singlet 770 , with its mass squared Mq remaining massive even in 
the chiral limit and behaving as 1 /Nq for Nq —> 00 , one needs to introduce the 1/Nc ex¬ 
pansion together with conventional chiral expansion consisting of soft momentum squared p 2 
and light-quark masses m q , in order to strictly establish a consistent power-counting [38l[39] . 
The simultaneous triple expansions on the soft momentum squared, light-quark masses and 
1/Nc, i.e. 0(5) ~ 0(p 2 ) ~ 0(m q ) ~ 0(1/Nc), will change the hierarchy of some operators 
assigned by the conventional chiral power counting. Typical examples are L^u^u^ 1 ) (x+) an d 
Lr > (u il u^x+) from the light flavor SU(3) yPT. Although the operators accompanying L 4 and 
L 5 belong to the same order in the SU (3) case, their orders are changed in the triple-expansion 
scheme when the singlet 770 is introduced [38], since they have different orders in 1 / Nq ■ Com¬ 
pared to the SU (3) octet case, new operators in the U (3) case with arbitrary number of the 
building block X = (In det U) oc 77 q term can appear [39], with U the matrix of the pNGBs. 
However, from large Nq argument, one more X field introduces one more order of the 1/Nc 
suppression. Therefore, the new operators with extra terms of the X field are more suppressed 
by 1/Nc- On the other hand, from the practical point of view, the operators with the X term 
exclusively contribute to the processes involving rj and rf mesons, which are only relevant to 
the coupled channels of the and pNGBs scattering. While the present lattice simulations 
in the coupled channel case are still rare and bear large uncertainties, as shown in later dis¬ 
cussions, it is still not possible to make sensible determinations of the LECs accompanied by 
the new operators with the X term. Moreover, since the singlet 770 predominantly contributes 
to the physical rf state, the contributions from the new operators with an X term will mainly 
enter the channels involving the rf meson, which have thresholds that are obviously higher than 
in the other channels. This further indicates the irrelevance of the operators with an X term 
in the present work. So we shall stop including any operator beyond the NLO discussion in 
Refs. [ T2 [1 13 il 1 8ll 19 j . which also enables us to make clear comparisons with the previous results. 

The LO Lagrangian that describes the interactions between the charmed-meson triplet 
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( 1 ) 


V = (D°, D + , Df) and pNGBs takes the form 

d£\ = v^vv^ - m 2 d vv +, 

which also provides the kinetic terms for the charmed mesons. The quantity Md stands for 
the mass of the charmed mesons in the chiral limit. The covariant derivative T>^ acting on the 
charmed mesons V is given by 


where 


V li V = V(S li + T] t ) , 


+ r M )pt, 


( 2 ) 
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In the above equations, r^ and denote the right- and left-hand external sources. F denotes 
the weak decay constant of the pNGBs in the chiral and large Nc limits, with the normalization 
F n = 92.2 MeV. A clear difference in the present discussion, compared to the previous ones in 
Refs. [T2l - H6lfT8l[l9l[2T] . is the inclusion of the singlet r ]o in pNGB matrix given in Eq. (|3|h 
The NLO Lagrangian that generalizes Eq. m reads [T2] 


^v\ = V (“Mx+) - ^iX+ + h 2 (u- h^u^) V ] 

+V^V (h 4 (uX) - h 5 {u», «"}) V v V t, (4) 


with 


X+ = u t xu t + ux^u , 11 ^ = i{v) (d^ - ir^u - u(d M - ui^u 1 } , (5) 

where \ = 2 B(s + ip) includes the scalar (s) and pseudoscalar (p) external sources. The 
quantity B is related to the light-quark condensate via (0|ifq^lO) = — F 2 B5^ at leading order. 
The light-quark masses in yPT are introduced by taking (s + ip) = diag(m, rh, m s ), with rh the 
average up and down quark mass and m s the strange quark mass. Isospin violating effects will 
not be considered throughout this work. Though in the strict triple expansion of U( 3) %PT 
the LECs in Eq. ([4]) belong to different orders, we simply quote them as the NLO ones as in 
Refs. [12L121] . since no other new operators with additional X held will be considered and in 
this way it enables us to make clear comparisons with the results in the literature. 

The relevant chiral Lagrangian with only the pNGB d.o.f is given by j33H35li38ll39] 

C x = F {u ^) + E (x+) + F-MqX 2 , (6) 

with X = ln(detC/). At leading order, the physical r/ and rj' can be obtained through the 
diagonalization of the r/o and rjs states in Eq. (0 , and they are related by 


V 




Vs 

Vo 



(7) 
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with sg = sin#, eg = cos # and # the mixing angle. The masses of the p and the rf and their 
mixing angle 9 can be expressed in terms of the parameters in Eq. ([6]) 


M 2 

AV± rj 

M 2 


sin# 


Ml , r2 

~Y + m k - 


Ml 


+ m 2 k + 


VMo 4 - 

4M 0 2 A 2 

+ 4A 4 

3 

2 

V M 0- 

4M 2 A 2 

+ 4A 4 

3 
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1+ 32A 4 


( 8 ) 

(9) 

( 10 ) 


with A 2 = Mj(- —M 2 . We mention that it is enough for us to consider the LO chiral Lagrangian 
for the pNGBs in Eq. ©. The LO square masses of the pion and the kaon are linearly dependent 
on the quark masses via 


M 2 = 2 Brh , Mk = B(m + m s ). 


( 11 ) 


With the previous setup, it is straightforward to calculate the charmed mesons (D and 
D s ) and pNGBs ( 7 and rf) scattering amplitudes, which can be categorized into seven 
different stes of quantum numbers characterized by strangeness ( S ) and isospin (/). The 
general expressions for the scattering amplitudes D\{p{) + (j>i(p 2 ) D 2 {p 3 ) + ^> 2 (^ 4 ) with 

definite strangeness and isospin can be written as 


V, 


( S,I ) 




Clo 


(s — u) — 4Cq/io + 2C\hi — 2024 ^ 24 ( 5 , t, u) + t, u ) 


where s,t,u are the standard Mandelstam variables, and H 2 i{s 1 t,u) and H^(s,t,u) are 


( 12 ) 


H 2 4 {s,t,u ) = 2h 2 p2 ■ PA + ■ P 2 P 3 ■ Pa + Pi ■ pm ■ P 3 ), (13) 

H 35 (s,t,u ) = h 3 p 2 -pa + h 5 (pi -p 2 p 3 -Pa + Pi ■ PaP2 ■ P 3 ) ■ (14) 

The coefficients Cj can be found in Table [l] and some explicit expressions of Ci are relegated 
to the Appendix. 

Notice that we have additional channels for the (S,I) = (1,0) and (0,1/2) cases compared 
to the previous SU( 3) calculations [41 fl2lfT3lfl9] and the scattering amplitudes involving the rj 
meson are also different from those references. Nevertheless, if we set the mixing angle 9 to 
zero and use the Gell-Mann-Okubo mass relation M 2 = (4 M^ — M 2 )/3, we explicitly verify 
that our results in Table Q] confirm the expressions in Refs. [UH2]. While comparing with the 
formulas in Ref. m, our results do not agree with the item at the intersection of the column 
labeled by Ci and the row labeled by D s r] —>• D s r]. We also find that there are global sign 
differences between ours and those in Ref. m for the Dir —>• Drj and Dir —>• D S K amplitudes 
in the case with ( S,I ) = (0,1/2). 

To continue the discussion, we perform the partial wave projections of the full amplitudes 
in Eq. (11211 and the projection formula with definite angular momentum J is given by 

V( lDm^D 2 <j> 2 ( s ) = \\ x d COS P ^( cos <P) V D%l^D 2( j>2 ( S ’ COS ^)) > ( 15 ) 
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Table 1: The coefficients in the scattering amplitudes ^ 1 0 1 ; _ k d 2 0 2 ( s j*>' u ) i n Eq. (fl2l) . The 
channels are labeled by strangeness (S) and isospin (I). For the coefficients not shown explicitly 
in this table, we give their expressions in the Appendix. 
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where cp stands for the scattering angle between the incoming and outgoing particles in the 
center-of-mass frame, and the Mandelstam variable t is given by 

t(s, cos ip) = M 2 Dx + Ml 2 — — (s + M 2 Dl — Ml) (s + — M^ 2 ) 

_£ SV a(s ' ■ M l > A < S - M k • M V> ■ (w) 

with A (a, b , c) = a 2 +b 2 -\-c 2 —2ab—2bc—2ac the Kallen function. We will be only interested in the 
S-wave projection in this work, i.e. the J = 0 case. Therefore, the subscript J labeling different 
partial wave amplitudes Vj Di<j>i-*D 2 <t >2 ( s ) dropped in later discussions for simplicity. 

The perturbative scattering amplitudes at any finite order alone can not generate resonances 
or bound states. Unitarity must be taken into account to study the resonances or bound 
states in the two-body scattering processes. We use the unitarization approach that has been 
extensively employed to discuss the T>* 0 (2317) previously in Refs. P^fl3lfl6lfl8l lT9]. This 
unitarization method was also used to study many other important phenomenons in hadron 
physics, such as the light-flavor meson resonances |22|28:,;4Q], the A(1405) [41] . etc. The 
unitarized two-body scattering amplitude in this formalism is given by [TQ 

T(s)=[l-V(s)-g(s)]- 1 -V(s), (17) 


where V(s) stands for the partial wave amplitude in Eq. (U5|) and for simplicity the super- and 
subscripts have been dropped. The function g(s) collects the unitarity cuts generated by the 
intermediate two-particle states in question. The loop function g(s) is 


g(s) = i 


d 4 q _ 1 _ 

(27t) 4 (q 2 - Ml + ie)[(P - q ) 2 - M| + ie] 


= P 2 


(18) 


which can be calculated in a once-subtracted dispersion relation or in dimensional regularization 
by replacing the divergence by a constant. The explicit form of g(s) reads j40] 


ff(s) = 


1 


167T 2 
a 
2s 


a(p) + In 


M 2 D s — Ml + M 2 _ M 2 


L D 


_ 

Ml 


^ 2s 

2s [ ~ + ° r ) — hi (—s + M| — M'l + a) 

+ ln(s + M^ — Ml + <j) — ln(— s — M^ + M'l + cr)] ^ , 


(19) 


with a = l\{s,Ml,M 2 ) and p the regularization scale. The function g(s) is independent 

of the scale p and the explicit ^-dependence in Eq. (fl9l) is compensated by the /x-dependent 
subtraction constant a(p). To be specific, we set p = 1 GeV in the numerical discussions, in 
order to make a clear comparison with the previous works mm- In this way, the value of the 
subtraction constant a(p) determined in our fits correspond to its value at the scale 1 GeV. 

The formalism for the unitarized partial-wave amplitude in Eq. (1171) can be easily general¬ 
ized to the coupled-channel case, where the quantities V(s) and g(s) should be understood as 
matrices. The matrix elements for V(s) can be calculated using Eq. (|15L Further, g(s) in the 
coupled-channel case becomes a diagonal matrix, with its matrix elements given by Eq. ciD 
with the pertinent values for Md and M^. 


7 











3 Fits to lattice simulation data 

Up to now, there are no experimental measurements on the charmed meson and pNGB 
scattering. On the other hand, lattice QCD simulations are much advanced in this research 
field gH7j and can provide us with valuable information to constrain the interactions between 
the charmed mesons and the pNGBs. The useful information that lattice simulations provides 
are the Dcfr scattering lengths. In the channel with definite strangeness and isospin, the S-wave 
scattering lengths are related to the unitarized scattering amplitudes T(s ) in Eq. (fl7]) via 

= ~ g+ Mff,) ( Sthr )' D< /’- >D 0’ Sthr = (Md + M,/,) . (20) 

In order to make a comparison with the exsiting lattice-QCD data, the scattering lengths 
should be extrapolated to the unphysical quark masses via 


Mr = \Jm 2 k + 1 -M* , 

M d = \Jm 2 d + 2(2/10 + /*i)M2 , 
M Ds = , 


( 21 ) 


o o o 

which are obtained from Eqs. (H|), Q and ©. Here, Mr, Md and Md s denote the masses in 
the limit of M% (oc m) —> 0, but with the strange quark mass m s ^ 0 (the so-called two-flavor 

o o o 

chiral limit). We mention that the quantities Mr, Md and Md s do not represent their values 
in the (three-flavor) chiral limit, since they are still dependent on m s . Their forms read 


Mr = Bm s , 

Md = M D + 4/io Mr , 

Md s = M D + 4(/io + hi) Mr , (22) 

where Md is the chiral limit mass of the charmed mesons, see Eq. ©• Combining Eqs. m, 
m and ( 1221 ) , one can easily perform the chiral extrapolations by varying the light quark mass 
rh and the strange quark mass m s . The corresponding light-quark mass dependences for M v , 
M^t and the mixing angle 9 are obtained by combinig Eqs. (HU), dED, ® and (1101) . Compared 
to the previous works in Refs. |12U18lfT9] . we do not further make the expansion of M n inside 
the square roots in the right-hand-sides of Eq. (I21D in order to obtain Mr, Md and Md s ■ It is 
worthy noting that the present lattice simulations are performed with fixed strange- and charm- 
quark masses, while only the up/down-quark masses are varied. In addition, the fixed strange 
and charm masses are usually set at (slightly) unphysical values in lattice QCD simulations 
and, furthermore, different configurations may adopt different values, which finally lead to 

o o o 

different values for Mr, Md and Md s ■ In our case, the lattice QCD data are taken from two 

o o o 

collaborations, e.g. Ref. [2] and Refs. and the values of Mr, Md and Md s corresponding 

to different lattice configurations, as well as their values extracted from the physical case, are 
listed in Table [2] for easy comparison. 
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Ref. gj 

Ensemble (1) [5H7] 

Ensemble (2) |5HZ] 

Physical 

M k [MeV] 

560.8 

519.0 

482.8 

486.3 

M d [MeV] 

1940 

1538 

1631 

1862 

M Ds [MeV] 

2058 

1655 

1731 

1967 


Table 2: The masses defined in the limit M n —>• 0 in Eq. (1221) for the chiral extrapolations. 


Next we present our fit procedure. The unknown parameters are F , Mo, hi= 0 , 1 , 2 , 3 , 4,5 and 
the subtraction constants a. Since we only work at leading order in the pNGB sector Q, it is 
justified to approximate the value of F by the physical pion decay constant F n = 92.2 MeV. 
For Mo, we adopt the LO value Mo = 835.7 MeV determined in Ref. {32], which has taken 
into account the recent lattice simulation data on the r/ and rf masses [33j. The LECs ho and 
h\ are determined by the masses of the D and D s mesons via Eq (f2ll) and ho by the slopes of 
lattice QCD data for the masses of the D and the D s while h\ is fixed from the physical mass 
splitting between D and D s . The resulting numerical values can be found in the first two rows 
in Table [3j 

Further, redefinitions of the remaining LECs, i.e. /i«= 2 , 3 , 4 , 5 j are introduced in order to 
reduce the correlations during the fitting process, as done in Refs. sum. These take the form 

/124 = h 2 + h'i , /135 = ho + 2 /ig , h^ = h^Mp , h'§ = h^Mp , (23) 

with My) = (M|) hyb + M]/ liys )/2 the averaged physical masses of D and D s . The redefined LECs 
are dimensionless and used as fitting variables. A common subtraction constant a, occurring 
in the loop function in Eq. (1191) , is adopted for all the possible channels @. Then the redefined 
LECs hi in Eq. (1231) and the subtraction constant a are fitted to the lattice scattering lengths 
provided in Refs. @H3. We stress that our fits are done using directly the lattice QCD data for 
the various masses, rather than those derived from Eq. (|21|) together with the masses shown 
in Table (2j We mention that our treatment of the masses in the fits is the same as that in 
Ref. pE] , but different from the one in Ref. [12], which used the masses resulting from the chiral 
extrapolations in Eq. (1211) . For reference, the correlation coefficients between the fit parameters 
for Fit-6C from the CERN MINUIT minimization package are shown in Table [H It is worthy 
noting that the correlations of the different lattice simulation data are not considered in our 
fits and hence the correlations of the fit parameters might be underestimated. 

Most of the data from the present lattice simulations are for the single-channel cases, e.g. 
(S,I) = (—1, 0), (—1,1), (2,1/2), (0, 3/2) and the D s n channel in the ( S,I ) = (1,1) case [3]. 
The lattice simulations for the ( S, I) = (1, 0) case, where the D* 0 (2317) appears, are still quite 
limited. Only two data points are given in Refs. [5H7]: one is determined with a Nj = 2 
simulation and the other is obtained with Nj = 2 + 1. To evaluate the effect of these two 
data points from the ( S,I ) = (1,0) simulations in the fit, we perform two different types 
of fits. The first one is a five-channel fit, denoted as Fit-5C for short, where the data for 

1 In general, the subtraction constant a(/r) should be channel dependent. However, if distinct subtraction 
constants are adopted in different channels, our fits tend to be unstable. This indicates that the present lattice 
simulation data are still not precise enough to discriminate the individual subtraction constants in different 
channels. 
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LEC 

Fit-6C 

Fit-5C 

Table V 0 

ho 

0.033* 

0.033* 

0.014* 

hi 

0.43* 

0.43* 

0.42* 

ft -24 

—0 ri + 0 - 05 
u ' ia -0.06 

—0 13+ 0 - 05 
u ' lo -0.06 

—0 in + 0 - 05 
u ' iu -0.06 

h'35 

oq +0.06 

u * zo -0.06 

0 24 +0 ' 12 
u ' zq: -0.12 

0 25 -1-0 ' 13 
u ' zo -0.13 

K 

—0 21 +0,29 
u - zl -0.27 

—0 IQ - * -0 - 32 
u ' iy -0.31 

fl qq +0.35 
U * OZ -0.34 

K 

A -'°- 0.19 

1 qq +0.57 
-*-•^-0.56 


o( A = 1 GeV) 

1 qc +0-06 

1 qc +0.07 
J -° O -0.09 

1 *°°-0.09 

X 2 /d.o.f 

12.27 -| r ) 

16-5 “ iiZ 

H 2 l - 1.23 

1.06 


Table 3: Fitting results of the parameters. The fitting results labeled by Table V are taken 
from Ref. [4] for comparison. Fit-6C and Fit-5C denote six- and five-channel fits, respectively. 
ho and h\ are determined by the masses of the D and the D s , not by the scattering lengths. 
See the text for details. 


Global 

^24 

^■35 

K 

K 

a 

h24 

0.983 

1.000 

-0.493 

-0.697 

0.135 

0.677 

h-35 

0.889 

-0.493 

1.000 

-0.165 

-0.546 

-0.761 

K 

0.971 

-0.697 

-0.165 

1.000 

0.501 

-0.001 

K 

0.954 

0.135 

-0.546 

0.501 

1.000 

0.780 

a 

0.981 

0.677 

-0.761 

-0.001 

0.780 

1.000 


Table 4: The correlation coefficients between different parameters result from the MINUIT 
program. 


10 
















the coupled DK(I = 0) channel are not included in the fit, e.g only the data points in the 
(S,I) = (—1,0), (—1,1), (2,1/2), (0, 3/2), (1,1) cases from Ref. [3] are considered. Compared 
to the five-channel fit shown in Table V from Ref. [3|> the coupled-channel effect for D s n is 
included here. To be more specific, we include the lattice QCD data corresponding to pion 
masses, 301 MeV, 364 MeV and 511 MeV from Ref. |4j in the Fit-5C case. The points for 
M 7l r > 600 MeV are excluded in the fits, since such a large value may challenge a reliable chiral 
extrapolation. The fit results for the parameters are given in the column labeled as Fit-5C in 
Table [3] and the reproduction of the lattice simulation data can be seen in Fig. |T] In the last 
panel of this figure, we make a prediction for the DK scattering length in the (S',/) = (1,0) 
case, where in order to make a clear comparison of the lattice data from Ref. [7] we have used 
the masses of the D and D s from that reference to make the plots, e.g. the results in the 
Ensemble (2) column in Table [2j Notice that we do not explicitly show the error bands of 
the predictions in the last panel of Fig. [H since the error bands are so huge that they almost 
cover the whole region of that panel. It is quite clear that our central predictions prefer the 
Nf = 2 + 1 simulation result over the Nf = 2 one. This inspires us to perform another type 
of fit, where we explicitly include the Nf = 2 + 1 datum from Refs. 0E] in the fit. This fit 
will be named as Fit-6C in later discussions, as in this case we have six channels in the fit. 
The fit results of the LECs and the subtraction constant are compiled in the column labeled 
as Fit-6C in Tableland the comparisons between our chiral extrapolations using Eq. (1211) and 
lattice QCD data are displayed in Fig. [2j Due to the inclusion of the coupled channel effects, 
improvements are achieved in the sense that the error bars become smaller, especially for the 
results in Fit-6C. Nevertheless, our results and those in Ref. [3] are compatible within the 
uncertainties and the resulting values for LECs are quite similar, see Table [3] for comparisons. 
The corresponding descriptions of the lattice QCD data are all of good quality with a y 2 per 
d.o.f around 1.0, which can also be seen from the plots in Fig. |T]and Fig. [2j In Figs. |T] and [2], 
the shaded bands represent the variation of the scattering lengths with the LECs whose values 
vary within a 1 -a uncertainty. 

Our prediction for the chiral extrapolation of Dix scattering length with (5,/)=(0,1/2) is 
given in Fig. [3l where we have used the masses in the column labeled by ’Physical’ in Table El 
In addition, with the LECs obtained in Table [3] and masses specified at their corresponding 
physical values, we predict all the physical S- wave scattering lengths in Tabled 

4 Pole analyses on the charm scalar states 

4.1 Pole contents in the physical situation 

Experimentally observed bound states and resonances correspond to the poles of the partial- 
wave amplitudes in physical and unphysical Riemann sheets (RS), respectively. In this section, 
we utilize the previously determined LECs and subtraction constant to study the pole contents 
in the charmed mesons and pNGBs scattering amplitudes. In our formalism, different Riemann 
sheets are characterized by the sign of imaginary part of the loop function g(s) in Eq. (|19l) . 
Each g(s) is associated with two sheets. By default the expression in Eq. (1191) defines the 
first/physical Riemann sheet. To reverse the sign of the imaginary part of the g(s) function in 
Eq. (|19l) . one then analytically extrapolates to the unphysical Riemann sheet. For scattering 
processes in question, different Riemann sheets can be accessed by choosing the proper signs 
of the imaginary parts of the Qi(s) functions, with gi(s) calculated from the two intermediate 
states in the z th -channel. The first Riemann sheet will be labeled as (+, +, +,...). The second 
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M [GeVl M k [GeV] M [GeV] 


Figure 1: (Color online) 5-channel Fit (Fit-5C). The bands represent the variation of the 
scattering lengths with the LECs varying within l-er uncertainty. The black solid circles with 
error bars stand for the lattice QCD data from Ref. [3J. The open and solid squares with error 
bars denote the Nf = 2 and Nf = 2 + 1 lattice QCD data from Ref. [6], respectively. The blue 
solid line in the (S = 1,1 = 0) panel is our prediction, see the text for details. 


and third sheets can be accessed by reversing the signs of the imaginary parts of the g(s ) 
functions defined at the first and second thresholds, which will be denoted by (—, +, +,...) and 
in order. Besides the pole positions, the residues for a given pole calculated in 
the partial-wave amplitudes in different channels can also provide useful information, since the 
residues correspond to the coupling strengths between the pole and its interacting Dcj) modes. 

Since the results in Table[3]for Fit-5C and Fit-6C are quite similar, we only present the pole 
analyses for the Fit-6C case in the following. Nevertheless, we stress that we have explicitly 
verified that the poles and their residues from the two fit results look indeed very similar. 
The pole positions and their residues in different channels, which are calculated by using the 
physical masses for the charmed mesons and pNGBs, are summarized in TableEJ Some remarks 
about the comparisons of our results for the charm scalar meson spectra with the previous ones 
in literature and experimental measurements are in order. 

A virtual-state pole at around 2.3 GeV is found in the single channel with (S', I ) = (—1, 0). 
In the (S, I) = (0,3/2) channel a very broad resonance with the width around 500 MeV is 
obtained. Such a broad resonance is difficult to be verified in experiments. 

For the (S, I) = (1,1) case with two coupled channels, we find one broad 2nd-Riemann- 
sheet pole above the DK threshold and one broad 3rd-Riemann-sheet pole below the DK 
threshold. The former one is mainly coupled to DK and the latter is coupled more or less 
equally to both D s n and DK. Our results are different from those in Ref. m, especially that 
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Figure 2: (Color online) 6 -channel fit (Fit- 6 C). The bands represent the variation of the scat¬ 
tering lengths with the LECs varying within 1-a uncertainty. See Fig. Q] for further details. 


our width for the 2nd-Riemann-sheet pole is much larger than that in the previous reference. 
These differences are presumably caused by the fact that the results in Ref. m are obtained 
replying on the preliminary scattering length data from Ref [33], which somewhat differ from 
the final published data that we used in this work, given in Ref. [3], It is quite interesting to 
point out that an enhancement in the D°K + invariant mass around 2350-2500 MeV was very 
recently observed [35], which may be an evidence of the existence of the pole with ^/s po i e = 
2466^27 - *271 ±5 in Tabled 

In the (5,/) = (1,0) case, a bound-state pole at 23211® is found and it corresponds to the 
physical Z)* 0 (2317) state. We further verify that the D* 0 ( 2317) is most strongly coupled to the 
DK system, as can be seen by its residues in Table [ 6 ] In addition, we also find a virtual-state 
pole just below the DK threshold on the second Riemann sheet. It is interesting to stress that 
the virtual pole is even closer to the DK threshold than the bound Z?* 0 (2317) state. 

In the (S,I) = (0,1/2) channel, there are two poles: ^/s po i e = 2114l| — illll® and 
v /s^I7 = 247311- i 1401®, located on the second and third Riemann sheets, respectively. The 
second-Riemann-sheet pole is consistent with the observations made in by Refs. jl4Hl6lfl9] . 
This pole can be regarded as a broad Dq (0, 1) state strongly coupled to the Dir threshold. On 
the other hand, the third-Riemann-sheet pole, which is most strongly coupled to D S K, has a 
much broader width compared to the previous results mmm and our determination of its 
mass and width is close to those of Dq(2460) announced by PDG [36] . 
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M. [GeV] ' M [GeV] 

Figure 3: (Color online) Predictions for the Dn —> Dtt scattering lengths in the (5 = 0, I = 
1/2) channel. The red lines with the cyan grid bands represent the chiral extrapolations of 
the scattering lengths with the masses specified by the last column in Table [2] and the LECs 
varying within their 1 -a uncertainties shown in Tabled) The open black square with error bar 
denotes the Nj = 2 lattice QCD data taken from Ref. [5]. The magenta solid square with error 
bar is our estimate based on the masses of Ensemble (2) in Table [2] and LECs in Table [3) 


(s, I) 

Channels 

Fit- 5 C 

Fit-6C 

(-1,0) 

DK ->• 

DK 

i 97+O.49 
- 0.36 

1 26 +0-46 

1 .ZU -032 

(-1,1) 

DK — > 

DK 

—0 21 +0 ‘ 02 
u - z± -0.01 

—0 21 +0 ' 01 
u - zl -o.oi 

(2,1) 

D S K -»■ 

D S K 

— 0 1Q +0 - 01 
u - ±y -o.oi 

—0 1Q +001 
u * iy -0.01 

(o,|) 

Dtt — > 

Dtt 

—0 101+ 0 003 
U.1U1_0 003 

—0 ioi+°- 001 
u.iui_ 0 .ooi 

(1,1) 

D s n — > 

D s tt 

0.004±g;ggi 

o.oo 4 ± 8 : 88 i 


DK -> 

DK 

o.o6±8;8i + * o.i7±8;8? 

o.o 6 ± 8;81 + < o.i7±8:8i 

(1,0) 

DK 

DI< 

-0.92ig;S 

- 0 . 89 t°;?o 


D s r] -*• 

D s r] 

-0.27±8;gl + * 0.03±g;gi 

-0.27±8:81 -H i 0.03±g;81 


D s rj' — > 

D s rf 

—0 22 +0 ' 03 -1- 7 0 ni + 0 - 01 
u.zz-o oi 1 u - ul -o.oi 

—0 22 +0 ' 01 4 - ? 0 Ol+ 0 - 01 
u - zz -0.0l ' ^ u - ui -0.0l 

(0,|) 

Dtt — > 

Dtt 

0 . 35 l°;°o 2 

0-3518:21 


Dr] 

Dr] 

0.02i° 0 ;°4 + *0.03±°; 0 0 ? 

o.o2i8;88 + i 0.0318:81 


D S K -4 

D S K 

-o.o 5 ± 8 : 8 S + i o. 35 ± 8;83 

-o.o 5 l 8;81 + * 0 . 3518:81 


Dr/ — >• 

Dr/ 

0-i6±8;it + * o.05±8;gi 

o.34l8:?i + *o.o4i8:S2 


Table 5: Predictions of the scattering lengths using the parameters from Table [3] together with 
the physical masses for the charmed and light pseudoscalar mesons. 
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(s, i) 

RS 

y/Spoie [MeV] 

(Residue] 1/2 [GeV] 

Ratios 

(-1,0) 

II 

23331“; 

7A5t 3 dt(DK) 



(0,1) 

II 

20331° - * 251I3 

6.64!°$(D7t) 



(1,1) 

II 

24661^ - * 27llg 

6-95!8;jj?(£U) 

1.72l°;{ 2 (DA:/D s 7r) 



III 

2225lg 2 - * 17811? 

7-35lg 

0.80t° o °ot(DK/D s 7T) 


(1,0) 

I 

23211® 

9.301 °o°A(DK) 

0.771 ° 0 f 2 (D sri /DK) 

0A3t° o il(D s r,'/DK) 


II 

235611 

2.851 ^(DK) 

0m^(D aV /DK) 

0.38! %£(D s rf/DK) 

(0,1) 

II 

2114 I 3 — i ml® 

9.66l°ii(i?7r) 

0.3l!°;°°(D ?? /D7r) 

0.49l° : °S( J D7//D7r) 

0.46l°;° 2 (D s A7Dvr) 


III 

2473I22 - * 1401? 

5.36!°!° (Dtt) 

1.09l°;°°(D ?? /D7r) 

1.12±%\l(Dr//Dir) 

2.12l° : °°(D s 77/Dvr) 


Table 6: Poles and their residues based on Fit-6C in Table [3j Physical masses for the charmed 
and light pseudoscalar mesons are used to obtain the results in this table. The Riemann sheets 
on which the poles are located are indicated in the second column. In the last column, we give 
the ratios of the residues with respect to the first threshold. 

4.2 Pole trajectories with varying pion mass 

Due to the limited experimental measurements of the charm scalar mesons, lattice simula¬ 
tions provide a possible way to verify the poles from our analyses. Therefore it is interesting 
to further probe the pole trajectories with varying pion mass M n , which can be useful for 
comparisons with future lattice results. From the theoretical point of view, the M n trajectories 
of the various poles can offer us further insights into the properties of different hadron states, 
as discussed e.g. in Ref. [23]. Before investigating the M n trajectories of the various poles, we 
first show how the masses of the K , rj, rj',D, D s , their thresholds and the mixing angle 6 change 
with varying pion mass in Fig. [4] We point out that in the discussions below we have used the 
physical strange quark mass, i.e. the values in the last column in Table [2j All of the masses 
considered here are increased when enlarging the pion masses, and the mixing angle 6 tends 
to increase as well. Around M n ~ 700 MeV, it is clear from Fig. [4] that the masses of D and 
D s become equal and the mixing angel 6 approaches zero. This phenomenon is not a surprise, 
since when M n ~ 690 MeV, the pion and kaon masses turn out to be degenerate, indicating 
the exact 517(3)-flavor symmetry in that region. 

The M n trajectories for the D* 0 (2317) meson can be seen in Fig. [5j In the left panel, 
together with the varying thresholds Mu + Mk, we give the results for the pole positions 
of D* 0 (2317) on the first Riemann sheet, which are identified as its masses. In the right 
panel, we show the binding energies which are the gaps between the thresholds and the pole 
positions. The red solid lines in Fig. [5] stand for our final results of the 17(3) chiral theory with 
realistic descriptions of i] and rf , while the blue dashed lines labeled as 517(3) are calculated 
by approximating our final expressions to mimic the SU (3) case, which are obtained by setting 
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Mq —>• ooH. In this limit, the heavy singlet rj 0 and the octet % will be decoupled, and one can 
identify the octet rjs as the physical 77 , as done in 577(3) yPT [37]. In addition, when Mq —>• 00 
the effects from the channels with the heavy singlet 7/0 will be much suppressed due to their 
far distance form the thresholds. The small variances between the U( 3) and SU( 3) lines in 
Fig. [5] indicate that the final results for D* 0 (2317) are not very sensitive to the treatment of 
the 77 and the rj' at least for low pion masses. This also explicitly justifies the previous study 
of D* 0 (2317) in Refs. [I2ffl6lfl8lfl9] based on the SU( 3) treatment of 77 . 

The most important conclusion obtained from Fig. [5] is that the D* 0 ( 2317) always stays a 
bound state below the DK threshold for a wide range of M n . In contrast, the M n trajectories 
of the pole around 2.1 GeV on the second Riemann sheet with (S,I) = (0, 1/2) in Table [6] are 
found to be quite complicated, as shown in Fig. [6l By increasing the value of M n , we first see 
that both the real and imaginary parts of this pole tend to decrease on the second Riemann 
sheet. At some point around M n ~ 2M| lhys , the real part of this pole becomes lower than the 
threshold of Dir, but its imaginary part is still nonzero. When M n equals to 288 MeV, the pair 
of resonance poles meets at the point below the Dir threshold on the real axis and becomes 
two virtual states on the second Riemann sheet. If we further increase the value of the pion 
mass, one of the virtual poles goes further away from the threshold towards to the left side 
on the real axis, while the other one moves closer to the threshold on the right side and it 
becomes a bound state on the first Riemann sheet when M n >336 MeV. If we keep increasing 
the value of M w , both the virtual and the bound state move further away from the threshold 
to the left side. It is very interesting to point out that the behavior for the broad charm scalar 
pole around 2.1 GeV with ( S,I ) = (0, 1 / 2 ) looks quite similar to the one of the a resonance 
/o(500) discussed in Refs. [ 221123] . 

The pole trajectories with varying M n for the charm scalar pole around 2.4 GeV in the 
third Riemann sheet with ( S,I ) = (0,1/2) are shown in Fig. [71 Compared to the trajectories 
in Fig. [6l the dependences of the pole around 2.4 GeV with are much weaker and only 
mild changes of the pole positions are observed over a wide range of pion masses from 100 to 
700 MeV. 

4.3 Nc trajectories for the charm scalar mesons 

The Nc trajectory of a resonance/bound-state pole can provide us useful information about 
the inner structure of the particle. Extensive studies on the Nq trajectories for light-flavor 
resonances, such as /o(500), /o(980), p(770), ai(1260) etc, have been carried out in litera¬ 
ture [241430] . As discussions on the Nc properties of the heavy-light mesons are still rare, one 
of the key goals of this work is to fill this gap. We point out that in order to make a reliable 
study of the Nc trajectories for a given state, not only the Nc scaling of the couplings of the 
scattering amplitudes, but also the Nc running of the masses of the intermediate particles in¬ 
volved in the scattering, need to be carefully considered. As mentioned in the Introduction, the 
singlet 770 , with its large mass due to the QCD U(1)a anomaly effect, should be properly taken 
into account when discussing the Nc dependences, as the U(1)a anomaly is 1 /Nc suppressed 
when Nc —>• 00 and 770 then becomes the ninth pNGB in the chiral and large Nc limits. 

For the Nc behaviors of the LECs in Eq. ((H) , their leading-order Nc scaling can be obtained 
by counting the number of traces in the corresponding operators. We take h 0 and hi as examples 

to illustrate this. The operator accompanied by hi has the same number of traces as the mass 

_2 

operator in Eq. ©, therefore hi has the same Nc scaling as the bare mass squared M D . The 
2 Strictly speaking, one would have to readjust the LECs, too. 
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Figure 4: Masses for the pseudo-Goldstones bosons and the charmed D mesons, the LO 77 - 77 ' 
mixing angle and the thresholds as functions of M n . 


leading scaling of ho is one more power of 1/Nq suppressed since there is one additional trace 
in this operator. From the large Nc point of view, the mass of any qq meson behaves as a 
constant @7][j8! and, therefore, the leading Nc scaling for Mo, M n and Mk is 0(1). As a 
result, the leading Nc scaling of h\ is also 0 ( 1 ), while ho should be counted as 0(1/Nc)- 
Similar rules can be applied to other operators in Eq. ((4]) and we summarize their leading Nc 
scaling as follows 

hi, / 13 , h 5 ~ 0(1), ho,h 2 ,h4~0(l/N c ). (24) 

The Nc scaling of the pNGB decay constant F in the chiral limit is 0(y/Nc) [36[[3Tj and the 
singlet 770 mass squared Mg is counted as 0(1/Nc) [33H36] , For the subtraction constant a in 
Eq. (1191) . a reasonable assignment for its leading Nc scaling is 0(1) as argued in Ref. [28] and 
we will also use this. With these assignments of the Nc scaling for different parameters, we 
can calculate the Nc running for other quantities. 

In Fig. [5] we show the Nc scaling of the masses of the pNGBs and charmed D mesons, 
the 77 - 77 ' mixing angle and the relevant thresholds. A striking phenomenon from the first panel 
in Fig. El is that the 7 / mass significantly decreases when Nc is increased and it tends to be 
equal to the pion mass in the large Nc limit. Contrary to this, when increasing the values of 
Nc, the mass of the rj decreases from the beginning, but it still has a relatively large mass 
around 700 MeV in the large Nc limit. These behaviors can be analytically understood from 
the leading order mixing formulas in Eqs. ©, © and (1101) . In the large Nc limit, one has 
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Figure 5: Properties of D s q( 2317) pole masses when varying M n : U(3) v.s. SU( 3). 


Mo —>• 0 and in this limit one can easily obtain from Eqs. (|8j), (f9|) and (1101) that 

M, n = M n , M r] ’ = \J2M^ - M%, arcsin 9 = ~]j^> (25) 

which perfectly explains the results in the two panels on the left side of Fig. 0 The results 
are consistent with the findings in Ref. [39j. The masses of the charmed mesons depend 
on the two LECs ho and h\ . as explicitly given in Eqs. (1211) and (1221) . The final results for the 
Nc running of Mjj and Me s are given in the top right panel in Fig. 0 which basically behave 
like constants when varying Nc- As a result, the Nc running of the thresholds as given in the 
bottom right panel is mainly caused by the running of M r] and M„/. An important finding is 
that the order of the thresholds M+ M v and Me + Mx is reversed and the former becomes 
lower than the latter for Nc > 7. 

The pole trajectories of the D* 0 (2317) with varying Nc are given in the left panel of Fig. [9] 
for the leading Nc scaling of hi- In order to highlight the relative position of the pole and the 
DK threshold, we have normalized the units of the real axis in terms of Mo+Mx- The physical 
D* 0 (2317) state is represented by the bound-state pole at Nc = 3 on the first Riemann sheet, 
i.e. the left most point in the real axis in Fig.[9l We also find a virtual-state pole below the DK 
threshold on the real axis on the second Riemann sheet at Nc = 3, which is even closer to the 
threshold than the bound-state pole. The explicit values of the pole positions at Nc = 3 can 
be seen in Table [6l When increasing the values of Nc, we observe that both the bound-state 
and virtual-state poles approach to threshold, and they finally meet at the threshold around 
Nc = 6. By further increasing the values of Nc, the bound-state and virtual-state poles 
move into the complex plane on the second Riemann sheet and become a pair of resonance 
poles. At Nc = 7 we see a kink structure in the trajectories and by increasing Nc afterwards 
the resonance poles move deeper and deeper into the complex plane, with increasing real and 
imaginary parts. The trajectory of D* 0 (2317) for large values of Nc is clearly different from the 
behavior of a standard quark-antiquark meson in QCD, which would fall down to the real axis 
with the mass behaving as a constant and the width scaling as 1/Nc |36j. In the left panel of 
Fig. M we show the Nc trajectories of the poles in the (S,I) = (0,1/2) case by considering 
the leading Nc scaling laws of hi in Eq. OH) . Compared to the D* 0 (2317) case in Fig. [9], the 
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Figure 6: Trajectories of the (S = 0,1 = 1/2) resonance at around 2.1 GeV with varying Al n . 
n is defined by M n = n M% hys 


trajectories of the two poles in the (S,I) = (0,1/2) channel show a much simpler behavior. 
For large values of Nc , both the real and imaginary parts of the two poles tend to increase 
and they keep moving into the complex plane instead of falling down to the real axis. This 
behavior look similar to that of the -D* 0 (2317) for large values of Nc, indicating that the two 
charm scalar resonances with (S,I) = (0,1/2) do not seem to be a standard quark-antiquark 
meson at large Nc- 

The previous Nc discussions are based on using the leading Nc scaling laws of the LECs 
hi given in Eq. (1241) . In order to check the stability of our conclusions on the Nc dependence 
of the poles, we follow Ref. [28] to include the subleading Nc scaling effects for the hi. The 
explicit formulae read 


hi(N c ) — hi(N c — 3) x < 1 + 


hi(N c = 3) - hi(N c = oo) / 3 


hi(N c = 3) \N 


-1 , * = 1,3,5 (26) 


and 


hi(N c ) — jy-hi(N c — 3) x 1 + 


hi(N c = 3) - h^ or (N c = oo) / 3 


hi(N c = 3) 


N r , 


-1 , * = 0,2,4, (27) 


with the normalized hf [or defined by hi Nor (N c = oo) = ^ hi(N c = oo). The fitted results of 
hi in Table [3] can be considered as their values at Nc = 3. The quantities hi(N c = oo) in the 
previous formulae stand for their values at large Nc, which can be estimated from tree-level 
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Figure 7: Trajectories of the (S = 0,1 = 1/2) resonance at around 2.4 GeV with varying M n . 
n is defined by M n = n Mh hys 


heavy resonance exchanges [50]. The resonance-saturation predictions of the hi from Ref. [50] 
are 

h 0 (N c = oo) = 0 , h 2 (1V C = oo) = 0 , h' A (N c = oo) = 0 , 

h\(N c = oo) = 0.42 , hs(N c = oo) = 2.23 , h' 5 (N c = oo) = —1.45, (28) 

which can be compared to our Fit-6C results 

h 0 = 0.033 , h 2 = -0.08^;^ , h! A = -0.21+°;^ , 

hi = 0.43 , h 3 = 3.79±g;ii , h' 5 = — 1.78±g;i| . (29) 

As we can see, the resonance-saturation determinations are quite consistent with the Fit-6C 
results. 

With the above preparations, we are at the point to study the influence of the subleading 

Nc scalings, specified by Eqs. (I26l27p . on the determinations of the pole trajectories with 

varying Nc- The corresponding results together with the ones by only considering the leading 
Nc scalings, are shown in Figs. iQlandfTOlfor the H* 0 (2317) and the poles in the ( S , I) = (0,1/2) 
channel, respectively. Qualitatively speaking, our conclusions based on the leading Nc scaling 
behaviors of hi are not changed when including the subleading Nc effects. Taking the i7* 0 (2317) 
trajectories for example, the only change is that the meeting point for the bound and virtual 
poles is at Nc = 8 in the case that the subleading Nc effects are included, while in the leading 
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Figure 8: Masses for the pNGBs and charmed D mesons, the LO rj-r)' mixing angle and the 
relevant thresholds as functions of Nc from 3 to 30 in steps of one unit. 

Nq case it happens at Nc = 6. But the trend of the pole trajectories when increasing the 
values of Nc are almost the same in both cases. Similar conclusions can be also made for the 
poles in the (S, /) = (0,1/2) channel, as can be seen from Fig. [TOl 

The 1 /Nc corrections of the subtraction constant a contributes another source of subleading 
Nc scaling effects. Since the subtraction constant is introduced through the unitarization 
procedure, it is rather difficult to directly estimate its leading Nc value, in contrast to the hi 
situation. We simply vary the values of the scale fx in Eq. (fTilL but keeping the subtraction 
constant a fixed at its fitted value, to roughly estimate the subleading Nc scaling effects of a. 
Several different values of the /j, from 0.9 GeV to 1.2 GeV are used to study the pole trajectories 
for the D* 0 (2317) and the ( S , I) = (0,1/2) channel. Since we do not observe qualitative changes 
by using different values of [i and they look similar as the cases by introducing the subleading 
Nc scaling effects in the hi , we shall not show explicitly these plots in order to avoid overloading 
the manuscript with too many figures. 

5 Conclusions 

In this work we have calculated the scattering of pNGBs (tt, K,r),r)') off charmed meson 
(D or D s ) and then unitarized the perturbative chiral amplitudes to investigate the charm 
scalar mesons. Recent lattice simulation data have been used in our study to constrain the 
free parameters in the fits. In addition to the prediction of the Dn scattering lengths with 
varying pion mass in the (S,I) = (0,1/2) channel, we also give predictions for all of the 
relevant scattering lengths with physical masses. Careful and extensive analyses of the charm 
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Figure 9: Nq trajectories of the D* 0 (2317) pole. The numbers in the plots indicate the values 
of Nc- Left: results with leading Nc scaling laws in Eq. (j2H> ; Right: results by including the 
subleading Nc scaling effects in Eqs. (1261) and (1271) . 


scalar poles and their residues in different (S, I) channels have been carried out, which are 
summarized in Table [6j A virtual pole is found to be responsible for the positive scattering 
length in the (S,I) = (—1,0) channel. Two poles are observed for the (5,/) = (1,1) case: 
the one on the second sheet may explain the enhancement recently reported in Ref. [45]. The 
pole positions for the H* 0 (2317) and the channel with (S', I) = (0,1/2) are compatible with 
previous determinations, with the exception that our determination for the (S,I) = (0,1/2) 
pole at around 2.4 GeV has larger width, which is close to the PDG value. Both M n and Nq 
trajectories for the D* 0 (2317) and the two poles in the (S, I) = (0,1/2) channel are studied in 
detail. The M n trajectories of the pole around 2.1 GeV quite resembles the ones of the /o(500) 
obtained in Refs. [221123] . The Nc trajectories of the D* 0 (2317) show that this physical bound 
state becomes a resonance for Nc > 6. For large values of Nc, the trajectories of the D* 0 (2317) 
and the poles with (S,I) = (0,1/2) do not tend to fall down to the real axis, indicating that 
they do not behave like the standard quark-antiquark mesons of QCD. 
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Figure 10: N c trajectories of the poles in the channel (0,1/2) with Nq varying from 3 to 30. 
The symbol cf) in the horizontal-axis label is n for the 2nd RS, and r/ for the 3rd RS. Left: 
results with leading Nq scaling laws in Eq. (|24l) : Right: results by including the subleading 
Nc scaling effects in Eqs. (1261) and (1271) . 


Appendix: Some of the coefficients in Table |T] 

Here, we give the form of the coefficients in Table [I] that were not explicitely stated there: 

■*1,0 DK^D s r) M^(5 cq + AV2s e ) + 3M^cq 


C 


C. 


1 2y/3 

1,0 DK^Ds-n _ Co + 2\/2 So 


35 


D s rj^D s rj 

°0 


^-y1,0 D s ri—>D a ri 


^1,0 D s ri—>D a ri 
°35 


V3 


- Ml) + AV2c e s e (M 2 K - M%) + s 2 g (2M* + Mp 
3 

2{Mt-2Ml){V2c e + s e ) 2 


2(\/2 eg + so)* 


(30) 

(31) 

(32) 

(33) 

(34) 


Ci 


1,0 DK^DsTj' 


c. 


c, 


C, 


1,0 D a r)—*D a rf 


c. 


M\{\\f2cQ - 5 so) + 2 >M^sq 


2\/3 


1,0 DK—tD s r)' 

35 — 


Sq — 2y/2 Cq 

V3 


1,0 DsO^Dgri 1 _ 

0 


1,0 D s ri—tD a ri' _ 


2{Ml - M 2 K )(V2c 2 g - cqsq - yggj) 

3 

2(2 M\ - M*)(V2c 2 e - cqsq - y/2s$ 
3 

-2(V2cj - cose - \[2sj) 


35 


(35) 

(36) 

(37) 

(38) 

(39) 
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( 40 ) 


/^l,0 D s rf ^ D s rj r 
°0 


£< 1,0 DsTj 1 — >DsT] f 


^1,0 Dsrf-^Dar)' 
U 35 


„ 0,5 Dr/^Dr 7 

0,| Dg—^Dg 

L 1 

0,i Dg^Dg 
°35 


sg(4A/% - i\4 2 ) + 4 % /2c ege (M^ - M 2 K ) + c 2 (2M 2 K + M(j?) 

3 

2(M 2 -2M 2 K )(^2sg-c e ) 2 

3 

2(y / 2s5> — eg) 2 
3 ’ 

c e(4Mjf - Mj) + 4V2c ege (M| - M^) + ^(2M|- + M 2 ) 

3 

—M 2 (y/2se — Cq) 2 

3 ’ 

(V2sq - Cg ) 2 

3 


„0,i D s K^Dg 
°1 

„0,1 D a K^Dg 
^35 


0,i Dg^Dg' 
°0 

^0,i Dg-^Dg' 

„0,i Dg-^Dg' 
°35 


cg(5Mx - 3 M 2 ) + Ay/2sgM^ 
2^6 

— (2\/2se + eg) 

vs ’ 

2(M^ - Ml){V2c 2 - cgsg - V2s 2 ) 
3 

M 2 (-V2c 2 g + cggg + V2sg) 

3 

(\/2c^ - Cgsg - V^Sg) 


(41) 

(42) 

(43) 

(44) 

(45) 

(46) 

(47) 

(48) 

(49) 

(50) 


r —y0, ^ D S K — >Dg* 
L '1 

0,i D s K^Dg' 
°35 


(5 Mft — 3M 2 )sg — 4\/2 M^cg 

2y/6 

2y/2 Cg — Sg 

Vi - ’ 


(51) 

(52) 


^-<1,0 Dg'—tDg' 
°0 


^1,0 Dg'—tDg' 
L 1 


^-<1,0 Dg'—tDg' 
°35 


gg(4M| - M^) + 4V2c eg 0(M2 - M|) + c 2 {2M 2 K + M^) 

3 

-M 2 (V2cg + gg) 2 
3 

(V2ce + sg) 2 
3 


(53) 

(54) 

(55) 
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